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ABSTRACT 

We investigate the growth of tidal instabilities in accretion discs in a binary star 
potential, using three dimensional numerical simulations. As expected from analytic 
work, the disc is prone to an eccentric instability provided that it is large enough to 
extend to the 3:1 resonance. The eccentric disc leads to positive superhumps in the 
light curve. It has been proposed that negative superhumps might arise from a tilted 
disc, but we find no evidence that the companion gravitational tilt instability can grow 
fast enough in a fluid disc to create a measurable inclination. The origin of negative 
superhumps in the light curves of cataclysmic variables remains a puzzle. 
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1 INTRODUCTION 

Classical superhumps are large amplitude variations in the 
optical light curves of close binary systems. The superhump 
period is typically a few per cent longer than the binary or- 
bital period. The most plausible explanation for these mod- 
ulations is that resonance of gas orbits in the outer disc with 
the tidal potential of the secondary star causes the disc to 
become eccentric (Whitehurst 1988; Hirose & Osaki 1990; 
Lubow 1991a). Due to a combination of the m = mode of 
the tidal potential, and pressure and viscous effects, the ec- 
centric disc precesses slowly in the prograde direction. Light 
from the disc is then modulated on the period of the disc 
motion relative to the binary potential (Murray 1996, 1998). 

In this paper we investigate the origin of a second, 
more puzzling, modulation in the light curves of some cat- 
aclysmic variables. These 'negative superhumps' have a pe- 
riod a few per cent less than the binary period. Their origin 
is presently unknown, but Patterson et al. (1993) proposed 
that this was the signature of the precessional motion of a 
tilted disc. This is an attractive explanation of the observa- 
tions, but if correct it poses the difficult question of how such 
a tilted disc can arise in cataclysmic variable (CV) systems. 

In this paper we investigate various methods of exciting 
a tilt or warp in a CV disc. The next section is a summary 
of the observations of the new periodicity. In section 3 we 
discuss the conditions under which tidal instability might in- 
duce a disc tilt, which we investigate numerically in sections 
4 and 5. Section 6 discusses the differences between CVs 
and X-ray binaries such as Her X-l, where the existence of 



warped discs is well established. Section 7 summarises our 
results, and discusses what observations are required to fur- 
ther investigate the possible existence of tilted discs in CVs. 



2 OBSERVATIONS OF NEGATIVE 
SUPERHUMPS 

An optical modulation with a period slightly shorter than 
the orbital period, P or t>, has been seen in several cataclysmic 
variables (Patterson et al. 1993; Patterson et al. 1997). Here 
we briefly describe a few systems that exhibit these so-called 
'negative superhumps'. 

The nova V603 Aquilae has a high mass transfer rate 
which implies a hot accretion disc (Patterson & Rich- 
man 1991). These systems are likened to dwarf novae in 
permanent outburst, and by analogy with those systems 
are expected to have a high (~ 0.1) a viscosity parame- 
ter (Shakura & Sunyaev 1973). V603 Aql's orbital period 
(Porb = 199.0 minutes) is longward of the period gap but 
still short enough to suggest that the disc could maintain 
contact with the 3 : 1 eccentric Lindblad resonance and be- 
come significantly eccentric. Indeed V603 Aql consistently 
exhibits a photometric signal with period ~ 210 minutes 
that is thought to be the signature of an eccentric, precess- 
ing disc (Patterson & Richman 1991; Patterson et al. 1997). 
Hence V603 Aql is sometimes referred to as a permanent 
superhumper. 

Patterson et al. (1997) reported the appearance of a 
second signal, with a period = 193.0 minutes, in the 
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V603 Aql light curve. This 'negative superhump' coexisted 
with and became brighter than its better known counter- 
part. They found the positive and negative superhump pe- 
riods varied on similar time scales and were negatively cor- 
related. In other words the negative superhump period was 
shortest when the positive superhump was longest. Patter- 
son et al. (1997) hypothesised that the accretion disc was 
simultaneously eccentric, and tilted. The prograde preces- 
sion of the disc's semi-major axis (as viewed in the inertial 
frame) gave rise to the positive superhump signal, whilst the 
retrograde precession of the disc's line of nodes was respon- 
sible for the shorter, negative superhump periodicity in the 
light curve. 

The SU UMa type dwarf novae V503 Cyg (Harvey et 
al. 1995) and VI 159 Orionis (Patterson et al. 1995) have 
both displayed large amplitude negative superhumps during 
quiescence and normal outburst. In the case of V503 Cyg 
the normal, positive superhumps were only seen during su- 
permaximum, whereas in V1159 Orionis they persisted well 
beyond superout burst and into the next normal outburst. 
Negative superhumps are not limited to short period systems 
however. For example TV Col, with P or b = 5.5 hours, dis- 
plays a photometric period — 5.2 hours (Hellier 1993). 



3 THE TIDAL INCLINATION INSTABILITY 

A fluid disc in a binary potential is subject to both eccentric 
and tilt instabilites at the 3:1 resonance (Lubow 1992a). 
Since it appears well-established that positive superhumps 
are the signature of an eccentric, precessing disc, we first 
consider whether the companion tilt instability can give rise 
to a negative superhump signal. In Lubow's analysis (1992a), 
the growth rates of these instabilities are (for the idealised 
case of a narrow gaseous ring at the resonance), 



2.08 Q olb q 2 ^ 
0.0398 Q olb q 2 



(1) 
(2) 



where W is the width of the resonance at radius r rcs , q is 
the mass ratio (M secon< j a ry/-Mprimary), and Q OI b is the binary 
angular velocity. To develop a tilt via this instability, we 
require that the growth time of the tilt T lnc ~ 1/Ainc must 
be much less than the drift time for material through the 
resonance, r a( jv = VK/|tv|. Making use of v r = — 3^/2r, and 
assuming an a prescription for the viscosity, v = ac 2 /Q.K, 
we obtain, 



2 2 (H\ 
Q « 225 g [T ] 



(3) 



where H = c s /^k is the disc scale height, c s is the sound 
speed, and fix is the Keplerian angular velocity. For q = 0.25 
and H/r — 1/20, the constraint on a from this analysis is 
a <g; 0.2. Given that typical estimates for a in high-state 
dwarf novae are > 0.1 (Cannizzo 1993), this immediately 
suggests that exciting a tilt via the instability is difficult. 
However it is also clear that these estimates are rather crude, 
and so in the next section we use numerical simulations to 
explore whether a gravitational tilt instability can be excited 
in the conditions of CV discs. 



4 NUMERICAL METHOD 



4.1 SPH code 



A smoothed particle hydrodynamics (SPH) code designed 
specifically for accretion disc problems has been described in 
detail in Murray (1996) and Murray (1998). Modifications 
have been made to extend the code to three dimensions, 
and to include adaptive resolution (i.e. to allow for variable 
smoothing length h). For disc simulations, the standard SPH 
linear artificial viscosity term (as described in Monaghan 
1992) can be modified slightly in order to improve its effec- 
tiveness as a source of shear viscosity in the SPH equations 
of motion. The modifications are described in some detail 
in Murray (1996) and Murray (1998). By taking the SPH 
equations to the continuum limit, it can be shown that this 
artificial viscosity term generates a shear viscosity 



v — k£ c B h, 



(4) 



where h is the SPH smoothing length, and ( is the (dimen- 
sionless) artificial viscosity parameter. For the cubic spline 
in three dimensions, n is analytically found to be 1/10. 

In the simulations shown in Murray (1996) and Murray 
(1998) the smoothing length had a unique fixed value. In the 
calculation described here we allow h to vary with both space 
and time. When evaluating the pressure- viscosity interaction 
between two particles, i and j, we choose to symmetrise 
with respect to the smoothing length and simply use h — 
\ (hi + hj) in place of h in the momentum equation. The 
density is evaluated using the form that appears as equation 
(6.3) in Monaghan (1992), namely 



Pi = W(n - r } , hj). 



We use 



-1/3 



(5) 



(6) 



to alter particle smoothing lengths at every time step. Here 
K is a constant set at the start of the calculation to ensure 
that each particle i has a reasonable number of neighbours, 
Ni. Neighbour numbers were monitored closely throughout 
the calculation and averaged near 75. The variable h imple- 
mentation of the code was tested against a fixed h version. 



4.2 Disc tilt and eccentricity 

We follow the evolution of tilt and eccentricity in the disc 
simulations by decomposing the particle distribution into 
Fourier components with argument (19 — m£l OI bt), where I 
and m are integers, and 9 is azimuth in the inertial frame. 
As in Lubow (1991b) we define the strength of the (l,m) 
mode 



Sl 7 rn(t) — ('S'cos^os^.m (t) + S cos , Bin ,i, m (t) 

H - 5'sin,sin,i,m 5sin, cos, l,m (tf)i. 



(7) 



The disc tilt is the (1, 0) Fourier component of the vertical 
displacement of the particles. The component modes of the 
vertical displacement are defined as follows: 



2fiorb 



7rN(l + 6i, )(l + S m ,o) 
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Figure 1. Grey scale column density maps of the disc at the times shown. The same density scale is used in all frames. The Roche lobe 
of the primary, the inner boundary at r w( j, and the approximate location of the tilt resonance are marked as solid lines. 



J^ JV _ 1 z p (t) sm(l6 p ) cos(mfl olh t') dt', (8) 



where 9 P is the angular position of particle p and N is the 
total number of particles. S is the Kronecker delta. We use a 
superscript z to distinguish modes of the vertical displace- 
ment from planar modes. 

As in Murray (1996) and Murray (1998), we use the 
interstellar separation d, the binary mass M, and the inverse 
of the orbital angular velocity as the length, mass and 

time scalings for our simulations. 



Table 1. Summary of parameter settings used in the tilt calcu- 
lation 

Parameter c s q ( limn »* wd At 

Units (if! or b d d ^orb 

Value 0.05 3/17 1.0 0.02 0.05 0.01 

5 RESULTS 

The tidal warping of accretion discs was studied by Larwood 
et al. (1996). They used three-dimensional SPH simulations 
to investigate the warping, precession and truncation of a 
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disc that initially had a non-zero inclination to the binary 
plane. In the case of CVs we must determine whether a 
tilt can be excited in a disc that is initially co-planar with 
the binary. In the simulation described here we investigate 
whether such a tilt can be generated by the inclination in- 
stability of Lubow (1992a). 



5.1 Initial conditions 

The initial condition for this simulation was an annulus com- 
posed of 21 concentric rings of particles spaced Ar = 0.01 d 
apart. Each ring consisted of 9 layers of particles centred 
on the midplane, spaced Az — 0.01 d apart. The disc then 
initially contained 29592 particles. Thereafter mass was 
added to the disc midplane at the circularisation radius 
r C irc = 0.1781 d at the the rate of one particle per time step 

At = o.oi n-^. 

We assumed the disc to be isothermal, with the sound 
speed c a = 0.05 d£l OI b- A binary mass ratio q = 3/17 was 
used. In previously described simulations (Murray 1998) 
with this value q, the disc rapidly became eccentric. We 
set the maximum value for particle smoothing lengths to 
be /imax = 0.02 d. We used an open inner boundary at 
TVd = 0.05 d. Particles ending a time step with r < r w a 
were considered to have been accreted by the white dwarf 
and were removed from the simulation. Particles that either 
returned to the secondary or were flung to very large radii 
(r > 0.9 d) were also removed from the calculation. The pa- 
rameter settings used in the simulation are summarised in 
table [jj. The calculation ran for 800 (approximately 127 
orbits of the binary). 



5.2 Disc structure 

Figures [lj through ^ summarise the structure of the model 
disc. Fig. ^] is a sequence of surface density maps showing 
the initial condition; the short time scale response; the later 
evolution; and the final state of the calculation. As expected 
for this mass ratio, the disc rapidly expanded to reach the 
Roche lobe of the primary, and to overlap the 3 : 1 eccentric 
resonance. Although there is considerable strength in the 
eccentric mode by the end of the simulation (see below), it 
is the (2, 2) tidal mode that is most apparent in each of the 
frames. 

Fig. ^ plots the azimuthally averaged vertical density 
profile for the annulus 0.348 d < r < 0.352 d, for comparison 
against the density profile of the corresponding isothermal 
atmosphere. With this particle number and choice of disc 
sound speed we obtain acceptable resolution of the disc's 
vertical structure over almost an order of magnitude in den- 
sity. Vertical velocities were subsonic everywhere in the disc, 
further indicating that the vertical resolution is satisfactory. 
For a cooler disc, with c s = 0.02 df2 rb we were unable to 
adequately resolve the vertical structure. 

The radial structure of the disc in the latter stages of the 
simulation is presented in Figures ^| and [l|. From Fig. |^ it is 
apparent that particle resolution is inadequate at small radii 
r < 0.15 d. Exterior to this radius, the effective Shakura- 
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Figure 2. Mean density (scaled to the midplane density) as a 
function of height above the midplane for particles in an an- 
nulus 0.348 d < r < 0.352 d (filled points) for the simulation 

at time t = 446.50 f2 _ * . At this time there are 34,406 parti- 
orb 

cles in the disc. The appropriate isothermal atmosphere (p(z) = 
p(Q) exp(—z 2 /2H 2 ) where H = rcs/v^ is the density scale height) 
is shown as a solid line. 



is approximately constant at a — 0.1. We expect this to be 
a reasonable value for alpha for nova-like systems such as 
V603 Aquilae, which have shown negative superhumps. It 
might be too high (and, the disc too hot) for dwarf novae 
in quiescence. In Fig. ^ the mean height of particles above 
the disc plane is shown as a function of r. Out to a radius 
r ~ 0.4 d, this profile is consistent with the r 3 / 2 scaling of an 
isothermal axisymmetric disc. At larger radii however, the 
influence of the companion and the non-axisymmetry of the 
disc breaks the r 3 ^ 2 scaling relation. 



Sunyaev a parameter, 

1 n 

° ~ 10 ^ c7' 



5.3 Mode strengths 

Fig. |E| shows the eccentric mode strength in the disc, defined 
in equation (Q), as a function of time. For times t < 550 
ft~ fb the eccentric mode grew exponentially, with a growth 
rate of A = 0.0044 ± 0.0001 Sl olh . When equation ^ is cor- 
rected to allow for a broad disc (see Lubow 1991a), we obtain 
the analytic estimate for the growth rate, A t ~ 0.02f2 m -b. 
The simulation eccentricity therefore grew significantly more 
slowly than analysis predicted, and also more slowly than in 
previous numerical calculations (Murray 1996, 1998). How- 
ever, the discs described in those papers were colder and 
more viscous than the disc detailed here. Consequently the 
tidal modes in those simulations were much weaker than in 
this calculation. Lubow (1992b) numerically showed that the 
non-resonant tidal response of the disc weakens the eccen- 
tric instability, an effect not accounted for in equation (Q). 
Hence the reduced eccentricity growth rate in this simula- 
tion is consistent with our understanding of the resonance. 
In Fig. plwe plot the strength of the disc tilt as defined 
(9) in equation (0) . The variations seen in this plot are of very 
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Figure 3. Estimated value of the Shakura-Sunyaev viscosity pa- 
rameter a as a function of radius for the simulation at time 
t = 446.50 ^~ rb - At this time the disc contained 34405 particles, 
and the mean number of neighbours was 77.5. 
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Figure 4. Average height of particles above the midplane a as 
a function of radius for the simulation at time t = 446.50 f2~ rb - 
The scale height for an isothermal disc with c a = 0.05 dQ OI ^ is 
plotted as a dashed line. 



small (<g 1 degree) amplitude, and are consistent with noise 
in the simulation. No significant tilt is generated. We ran 
a short comparison calculation, taking the disc output at 
t = 634.50 Sl~\ and adding a one degree inclination to it. 
Mass addition, and all other parameters were kept as before. 
The subsequent evolution of the tilt mode strength is shown 
in Fig. (Q). With the addition of the tilt, j is almost two 
orders of magnitude larger than at any stage of the original 
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Figure 5. Eccentric mode strength as a function of time. 
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Figure 6. Strength of the (1,0) mode of the vertical displace- 
ment as a function of time. For graphical convenience we have 
actually plotted 10 5 X . The strength of this mode was only 



calculated for t > 400 Q 



simulation. The mode strength however then decayed with 
a time scale ~ 400 fi~ rb , which agrees well with the simple 
estimate Ma BC /M ~ 350 fi~ rb . This is as expected if the 
continued addition of mass simply dilutes any initial disc 
inclination. 



6 COMPARISON WITH X-RAY BINARIES 

The prototypical example of a warped, precessing accretion 
disc is that of the X-ray binary Hercules X-l. Recently it has 
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Figure 7. Strength of the (1,0) mode of the vertical displace- 
ment as a function of time for disc with perturbation = 
Zi + 0.02 r; cos(0;) applied at time t = 634.50 



-l 

orb ' 



been suggested that this, and similar systems, develop warps 
as a result of the radiative warping instability analysed by 
Pringle (1996; see also Maloney, Begelman & Pringle 1996; 
Maloney & Begelman 1997) . Here, we briefly summarise why 
the same mechanism does not work for CVs, except in a few, 
exceptional, circumstances. 

The radiative warping instability arises because repro- 
cessed radiation from a central luminous source is re-emitted 
non-axisymmetrically from an optically thick disc possessing 
a warp. Growth occurs if the disc is larger than some critical 
radius, r C rit, given for an isothermal disc (Pringle 1997; for 
the general case see Maloney, Begelman & Nowak 1997) by, 



2^ 



(10) 



where r s = 2GM/c 2 is the Schwarzschild radius of the ac- 
cretor, e = L/Mc 2 is the radiative efficiency of the flow, and 
r\ is the ratio of the 'vertical' viscosity (that acting to reduce 
the relative inclination of disc annuli) to the shear viscosity 
in the plane (see Papaloizou & Pringle 1983 for a formal 
definition). For accretion powered cataclysmic variables, 

~ GM in-* nil 
e» j ~ 10 , (11) 

and the instability criterion, equation ([U)|), requires r cr it pro- 
hibitively large, of order 10 14 cm - much greater than the 
actual extent of the disc. Steady CV discs are thus generally 
stable against radiative warping. 

Evading this conclusion requires either a higher radia- 
tive efficiency or an appeal to non-steady-state conditions. 
The former is possible if there is steady nuclear burning of 
accreted material on the surface of the white dwarf. The ra- 
diative efficiency in this regime is as high as e ~ 7 x 10 , 
leading to an estimate for r cr i t of the order of 10 1 cm. 
With the substantial theoretical uncertainties, for example 
in 77, this could lead to warping provided the disc was large 



enough and optically thick. However the requirement for nu- 
clear burning limits these possibilities to supersoft X-ray 
sources (Southwell, Livio & Pringle 1997) rather than lower 
accretion rate CVs. 

The wide range of transient behaviour seen in CVs sug- 
gests a number of scenarios where non-steady-state warping 
effects might be important. For example, the outer cool parts 
of a dwarf nova disc might be susceptible to warping from 
the radiation emitted by hotter inner regions that were de- 
caying from outburst, or from radiation from a cooling white 
dwarf heated during the high state accretion event. Unfortu- 
nately, while these scenarios produce a disc that is formally 
unstable by equation ( |i"o| ) , the time scale for growth is gen- 
erally much greater than the duration of the transient state. 
The time scale for growth is (Pringle 1996), 



Trad 



127rEr 3 f2xc 



which for highly optimistic CV parameters is, 
1/2 



r rad ~ 5 x 10 
X (l0 10 cm) 



M 



0.6 M Q ) VlOgcm- 2 
3/2 / L 

2 x 10 33 ergs- 1 



(12) 



(13) 



We therefore conclude that there is insufficent luminosity 
to generate a finite warp in this manner. Moreover there is 
no sign that the systems with observed negative superhump 
signals are universally prone to dramatic transient changes 
in accretion rate, as would certainly be required for this 
mechanism to operate. 



7 DISCUSSION 

We have found that the tidal inclination instability is too 
weak to generate a significant tilt in an accretion disc in 
the high state. However, neither our time scale analysis, nor 
our simulations can rule out such instability occurring in 
short period systems with quiescent discs. In a disc of very 
weakly interacting particles, vertical instability as opposed 
to tilt instability was relatively easy to excite at the 3 : 1 res- 
onance. Thus, although it seems unlikely that the inclination 
instability discussed by Lubow (1992a) can account for the 
negative superhump phenomenon in currently observed sys- 
tems, the instability might be important in low state dwarf 
nova, particularly those such as WZ Sge systems where the 
quiescent a parameter is thought to be very low (though see 
also Warner, Livio & Tout 1997 for a contrary view). Al- 
though the contribution to the observed emission from the 
disc is small in these systems, they are likely to be the best 
places to investigate tidally-driven disc tilts. 

Although this work does not categorically rule out the 
possibility of tilted discs in cataclysmic variables, their ori- 
gin is somewhat doubtful. The radiative torques speculated 
to be responsible for the better established warps in X-ray 
binaries definitely do not generally operate in CVs, leav- 
ing no clear mechanism that is known to be effective in 
creating a tilted disc. We note that of those systems that 
do show negative superhumps, TV Col is a known inter- 
mediate polar (Hellier 1993), and V603 Aql is a suspected 
one (Schwarzenberg-Czerny, Udalski & Monier 1992). That 
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then implies that the dynamics of the inner disc are signif- 
icantly influenced by the magnetic field. The details of the 
role a magnetic field might play in creating a warped disc 
are not clear, however we note that Agapitou, Papaloizou & 
Terquem (f997) have shown that strongly magnetised thin 
discs are unstable to bending modes. Further investigation 
of the correlation between negative superhumps and weakly 
magnetic CVs would be very interesting. 
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